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The micro-structural organisation of a granular system is the one most important determinant of 
its macroscopic behaviour. Here we identify the fundamental factors that determine the statistics 
of such micro-structures, using numerical experiments to gain general understanding. The exper- 
iments consist of preparing and compacting isotropically two-dimensional granular assembUes of 
polydisperse frictional discs and analysis of the emergent statistical properties of quadrons - the 
basic structural elements of granular solids. The focus on quadrons is also because they are central 
to statistical mechanics of granular matter. The dependence of the structures and of the packing 
fraction on the inter-granuar friction and initial state are analysed and a number of significant 
results are found, (i) An analytical formula is derived for the mean quadron volume in terms of 
three macroscopic quantities: the mean coordination number, the packing fraction and the rattlers 
fraction. Using a recently-derived equipartition principles for such systems, this relates directly the 
mean coordination number, the rattler-free packing fraction and Edwards compactivity. (ii) We 
derive a unique, initial- state-independent, relation between the mean coordination number and the 
rattler-free packing fraction, showing that the common wisdom in the engineering community that 
such a relation depends on the initial state is misguided. The relation is supported numerically for 
a range of different systems, (iii) We collapse the quadron volume distributions from all systems 
onto one curve and verify that they all have an exponential tail, as predicted by the statistical 
mechanical approach - an exponential form that relates directly to a Boltzmann-like factor, (iv) 
The nature of the quadron volumes distribution is investigated by decomposition into conditional 
distributions of volumes given cell order and we find that each of these also collapses onto a single 
curve, (v) We find that the mean quadron volume decreases with increasing inter- granular friction 
coefficients, an effect that is prominent in high order cells. We argue that this phenomenon is due to 
an increased probability of stable irregularly-shaped cells and test this by a herewith developed free 
cell analytical model. We conclude that the micro-structural characteristics are governed mainly by 
the packing procedure, while effects of inter- granular friction and initial states are details that can 
be scaled away. 



I. INTRODUCTION 

This paper discusses the structural characteristics of 
random packings of planar (2D) granular solids. The sci- 
ence of random granular packing has a very long history 
and statistical methods have been introduced into this 
field in the 60's [H, Q. In 1989, Edwards and co-workers 
[sl, 13 proposed a novel statistical mechanical approach, 
based on the entropy of configurations that assemblies of 
grains can take. In this approach, energy is replaced by 
volume, the Hamiltonian is replaced by a volume func- 



* [tmatsu@kz.tsuku ba.ac.jp"] 
^ lrbbll@cam.ac.uki 



tion and the temperature, characterising ensemble fiuc- 
t nations, is replaced by compactivity. This latter quan- 
tity has been elusive to measure and quantify from bulk 
quantities. Much work followed, aiming both to use this 
formalism for derivation of macroscopic quantities, as in 
thermodynamics, and to test it experimentally and nu- 
merically [EUrj. 

Fundamental to this approach is an identification of 
basic volume elements, which can be used as 'quasi- 
particles' for the volume function. Blumenfeld and col- 
laborators have proposed that these quasi-particles are 
quadrons [8|-[iq| - structural elements that tessellate the 
granular space and can be used to quantify the local 
structure [1, Ell, [H. Ball and Blumenfeld [8] used the 
quadron description to construct a stress theory of iso- 
static granular materials. In contrast to other definitions 
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of volume elements, such as Volonoi-based tesserations 
Q, quadrons preserve the connectivity information, since 
their construction is based on the force-carrying inter- 
granular contacts. They also allow for an unambiguous 
tensorial description of every structural element. A fur- 
ther advantage of this description is that, at least in 2D, 
there are as many quadrons as there are degrees of free- 
dom in the statistical mechanical structural phase space 
S [iol, [HI . This makes it useful for treating these volumes 
as degrees of freedom in their own right. 

It has been shown that the quadron volume proba- 
bility density function (PDF), P{Vq), is best understood 
in terms of its conditional PDF's P{Vq \ e), where e is the 
order of the cell associated with the quadron (see Figure 
[1]). Based on geometrical considerations, and supported 
by numerical results, Frenkel et al. argued that the PDFs 
P{Vq I e) should be independent of inter-granular friction 
and of the packing process (i.e. the history). A detailed 
analysis of this issue and an understanding of the shape 
of these PDFs from first principles is yet to emerge. 

Here, we present a set of numerical experiments on as- 
semblies of polydisperse discs, isotropically packed by the 
same process to different packing fractions and mean 
coordination numbers z. We investigate in detail the 
geometrical characterisation of the quadrons and the de- 
pendences of the PDFs P{Vq) and P{Vq \ e) on (j) and z. 
We examine the decomposition proposed by Frenkel et 
al. [11], focusing on the significance of cell shapes, and 
we show that the universality suggested by them must be 
augmented with a mechanics-based consideration. This 
is supported by a free cell model, whose results are com- 
pared with the numerical experiments. We conclude with 
a discussion of the results and the new insight that they 
provide into structures of random granular packs in me- 
chanical equilibrium. 



II. PROBABILITY DENSITY FUNCTION OF 
QUADRON VOLUMES IN STATISTICAL 
MECHANICS 

As in thermodynamics, one can define the entropy, 5, 
of a micro-canonical ensemble of granular assemblies in 
static equilibrium as the logarithm of the number of ways, 
M, that a collection of N grains fit into a given volume 
V, 



degrees of freedom (SDF), {r^}, for any grain configura- 
tion. B is a product of delta functions, constraining the 
possible configurations to those specified in the ensemble 
statistics, e.g. to assemblies in mechanical equilibrium, 
with a specified mean coordination number z. From this 
micro-canonical description, one can readily pass to the 
canonical one [14]. Assuming a Gibbs ensemble [sl, 
the occurrence probability of a configuration of volume 
is 

where X = {dS/dV)~^ is a Lagrange multiplier, called 
compactivity, that characterizes the statistical fluctua- 
tions and is analogous to the temperature. Zy is a volu- 
metric partition function 

Z„ = ^e-^-/^^ /'e-'^(W)/^e(M)d^-/ M (4) 

i 

To use the statistical mechanical formalism, one needs 
a set of elementary entities, whose summed volumes gives 
the volume of the system. The quadrons, discussed 
originally in [8[ and adapted to statistical mechanics in 
S ilPj , are structural elements that fit this role nicely. 
They can be regarded as the 'quasi-particles' of this 
formalism [il, [io|- One also needs to identify the SDF. 
In terms of these, the volume function is Vq {{t]})- It 
has been shown [ol, [io| that Nsdf is equal to the number 
of quadrons in 2D, but is larger in 3D[15|. 

An explicit illustration of the quadron in 2D is shown 
in Figure [TJ First, one defines the centroids of grains g 
and of cells c as the mean position vectors of the con- 
tact points around them, respectively. Second, the con- 
tact points around every grain are connected to make 
polygons, whose edges are vectors, f^^, that circulate the 
grain in the clockwise direction. Third, one extends vec- 
tors R^^ from the centroids of grains g to the centroids of 
the cells c that surround them. A quadron is the quadri- 
lateral whose diagonals are the vectors f^^ and R^"^. Since 
every pair gc corresponds exactly to one quadron, we 
shall denote a quadron by the index q for convenience. 
The quadron's shape is quantified by the following struc- 
ture tensor 



S{V,N) = log M (1) 
where M can be written as 

M = ^S{V-Vi)Gi^ f S{V -W) e d^^^f{r]}{2) 

Here 6 is the delta function, Gi is the number of configu- 
rations (microstates) in mechanical equlibrium having a 
total volume Vi^ W is the volume function - a function 
that gives the volume in terms of Nsdf internal structural 



= (5) 

and its volume is 

V, = ^\7»xR'i\=^{r^,Rl-r1^Rl) (6) 

It has been shown that, in 2D, there are as many 
quadrons as there are degrees of freedom [ol, |lo|, [III , Nq = 
Nsdf - This makes it possible to use the quadron volumes 
to span the A^^-dimensional phase space, V. We define 
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FIG. 1. Illustration of quadrons 
(shaded) and cells (yellow) in 2D. 



computed using Newton's second law. We postulate a re- 
pelling harmonic interaction potential, characterized by 
normal and tangential spring constants, kn and /cg, ac- 
tivated upon contact and overlap between discs. We set 
ks/kn = 1/4 in this study. 

Our aim is to study polydisperse systems and we have 
chosen the distribution of the radii of the discs to be 
log-normal distribution due to its wide use in civil engi- 
neering and soil sciences [HI, 



P{D) 



1 
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exp 



{\nD-\nDof 
2^2 



(10) 



where we assign I^o = 1- and a = 0.2, which give for the 
mean grain size and mode, respectively, D = 1.02 and 
D^ode = 0.961 (see Figure [2|). 



the 'density of states' g (v^ d^^V = g {{Vq}) Uq dVq as 

the number of states whose volumes are between V and 
V -^dV. 

The volume partition function can be then written as 



Na 



g[{Vq})Y[dVq^ 



(7) 



where the integration is understood as taken only over 
the systems that satisfy the constraints in 0({?7}) and 

Let us now consider a model system - an ideal quadron 
gas (IQG) - where the quadrons are uncorrelated. The 
general partition function ([7j) becomes 



(8) 



where giiV) is a one-quadron density of states [10]. The 
propability density function (PDF) of finding a quadron 



of volume between Vq and Vq 



dVq is 



PiVq) 



-vjx 



(9) 



In the following, we investigate this PDF in our nu- 
merical experiments and relate it to the cell structures. 
The combination of the exponential term and the den- 
sity of states, whose exact form is unknown, makes exact 
calculations impossible without additional information. 
Nevertheless, it is possible to gain insight from such a 
study, as we show in the following. 



III. NUMERICAL EXPERIMENTS 

Our numerical experiments have been carried out us- 
ing the Discrete Element Method (DEM) ([ii,[l3). The 
method consists of using an incremental time marching 
scheme, wherein the motions of 2D discs (grains) are each 
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FIG. 2. Grain diameter probability density func- 
tion. 

The packing protocol of our systems is as follows. 
First, we construct three random packs in a double pe- 
riodic domain, engineered to be on the verge of jam- 
ming. The packs consist of about 21400 ± 1000 discs and 
are made at packing fractions (j) = 0.76,0.82 and 0.84. 
The small variation in the numbers of particles is due to 
the different densities of the three initial configurations. 
These configurations are then used as initial states for 
the packing procedure; loose initial state (LIS), interme- 
diated initial states (IIS), and dense initial state (DIS), 
respectively. 

Once an initial state is set, we assign all the parti- 
cles the same friction coefficient, /i, and apply to the 
system a slow isotropic stress ac by changing the peri- 
odic length on both sides. We limit the applied stress 
to a level that corresponds to an average overlap discs 
of S = (Jc/kn = 10~^. No gravitational force is applied 
and the compression continues until the fluctuations of 
both grain positions (per mean grain diameter) and inter- 
granular forces (per mean average contact force) are be- 
low very small thresholds - 10~^ and 10~^, respectively. 
This procedure is carried out for each initial state at five 
different values of the inter-granular friction coefficients: 
/i = 0.01,0.1,0.2,0.5 and 10, giving altogether 15 assem- 
blies. 
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Using this procedure, we have computed the packing 
fractions, the mean coordination numbers, and studied 
the structures in these systems. For the determination 
of ^, we disregard 'rattlers', i.e. grains with one or no 
force-carrying contact. The results are shown in Figure 
[3l The upper and lower bounds of the coordination num- 
ber Zmax = 4 and Zmin = 3 corrcspoud to the isostatic 
states for discs with friction coefficients of /i = and oo, 
respectively. The three systems with fi = 0.01 converge 
into states with (z, ^) = (3.941, 0.840), (3. 944, 0.842) and 
(3.943, 0.843), respectively, which are very close to the 
ideal frictionless jammed state. However, with increasing 
friction coefficient, the difference between the values of z 
and (f) in the final states of systems started from different 
initial states increases (see Figure [3]). 




0.79 0.80 0.81 0.82 0.83 0.84 
packing fraction ^ 

FIG. 3. Coordination number vs. packing fraction for 15 
systems, generated from three different intial states, LIS, 
IIS, and DIS. Note the convergence for /i ^ 0. 

Turning to the analysis of cell structures in these sys- 
tems, consider the three examples shown in Figures |4]-[6l 
constructed from LIS, with /i = 10 and 0.01 (|5|), and 
from DIS with = 10 (|6]). While traces of the initial con- 
dition can be detected in both the first two systems, for 
higher friction one clearly ends up with typically larger 
cells and a correspondingly lower z. Note that the pack- 
ing fraction in Figures |4] and [6] are very similar (see also 
Figure [3]) , but the cells in |4] are noticeably bigger typi- 
cally. Also note the considerable number of rattlers in 
the large cells, an issue that will be discussed in detail 
below. To quantify the structural differences, we next 
study the quadron volume distribution. 

IV. STRUCTURAL CHARACTERISTICS 
A. Quadron volume statistics 




FIG. 4. Example of an assembly with 
/i = 10, generated from LIS. 




FIG. 5. Example of an assembly with 
/i = 0.01, generated from LIS. 




FIG. 6. Example of an assembly with 
/X = 10, generated from DIS. 



all PDFs of quadron volumes are plotted in Figure [7] for 
all friction coefficients and for LIS and DIS. The corre- 
sponding means are plotted against the mean coordina- 
tion numbers, for all initial states, in Figure [H All the 
points fall nicely on one curve regardless of initial state. 

On close inspection, this very weak dependence, if any, 
on initial state can also be observed in Figure [3 This 
implies that the dependence of the relation between v 
and z on the initial state, seen in Figure [3l can be made 
to disappear under the right parameterization. A clue 
to such a parameterization is the fact that the quadron 
volumes are computed ignoring the rattlers. Following 
this idea, we can derive the relation between v and z as 
follows 



The grain size distribution is the same for all packs, 
making it convenient to normalize the quadron volumes 
by the mean grain volume Vg, v = Vq/Vg. The over- 
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FIG. 7. PDF of the quadron volumes, normalized 
by the mean grain volume, v = Vq/Vg for all the 10 
systems generated from LIS and DIS. 
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FIG. 9. Mean coordination number z vs. 
rattlers- free packing fraction (f)\ 
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FIG. 8. Mean quadron volume v vs. mean coordi- 
nation number z for all initial states. 



it is the small differences between the curves that gives 
rise to the observed differences in the conventional plots. 
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(j)' z (1 — rir)4>z 

where r]r = N^/N is the rattlers fraction and (j)' is the 
packing fraction after removing all the rattlers. This el- 
egant result suggests that there is indeed a unique rela- 
tionship between z and the packing fraction, but only if 
the latter is corrected for the rattler fraction. This is ver- 
ified by a direct plot of these two quantities, Figure [9l It 
is also important to note that this relation is independent 
of initial state as Figure [9] shows. 

Comparing Figures [9] and [3l we see that the rattler- free 
packing fractions are much smaller than the packing frac- 
tions commonly reported in the literature. The removal 
of rattlers also makes sense on mechanical grounds, since 
it does not affect the force transmission in our systems 
(nor in any system experiencing no body forces). 

We conclude that the seeming sensitivity of z{(j)) to the 
initial state (e.g. Figure [3]), commonly seen in the litera- 
ture, stems directly from the variation in r]r. Plotting rjr 
as a function of z in our systems (Figure [10]), we see that 



FIG. 10. The rattlers fraction r]r vs. the mean 
coordination number z. 

It is important to comment here that this result is not 
really universal. It is correct for the specific packing pro- 
tocol used here. While we believe that, for any given pro- 
tocol, the plots of z{(j)^) would collapse onto one curve, 
there is no reason that this curve should be universal. 
In other words, we expect different protocols to display 
different z{(j)') curves, probably depending on the rate of 
rattlers generation. However, these results do provide a 
universal protocol-independent insight - the relation be- 
tween the z and the packing fraction is directly linked 
to the mechanical stability of the structure and the way 
that forces are transmitted. We will discuss this insight 
more in the concluding section. 

Turning to consider the PDFs of the quadron volumes 
more closely, it is natural to expect the mean quadron 
volume to increase with /i, simply because the cells get 
bigger. To get insight into the shapes of the PDFs, we 
scale the quadron volumes by their means, u = Vq/Vq = 
v/v. This simple scaling is sufficient to collapse all the 
PDFs almost perfectly onto one curve for all the systems. 
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FIG. 11. The PDF of the normalized quadron vol- 
ume u = Vq/Vq for for all the 10 systems generated 
from LIS and DIS. 



independent of friction and initial state (Figure [TT]). 

Moreover, the collapsed PDF has an exponential tail, 
which is a signature of the Boltzmann-like factor of equa- 
tion [9l Our best fit to the collapsed curve has the follow- 
ing form 



P{u) = f{u)e 



(13) 



where a is a slope of the tail and f{u) is a rational func- 
tion. Comparing this expression with equation [9l we 
identify 



X = V,/a 



(14) 



where we used u = v/v = Vq/Vq. This implies that the 
compactivity X is only a function of Vq, in agreement 
with the equipartition principle derived recently in [l5|. 
Intriguingly, combining this observation (or the equipar- 
tition principle of [l5|) with equation ([12]), we find a ex- 
pression relating directly the mean coordination number 
the compactivity X and the rattler-free packing frac- 
tion 6' 



a 



^'ZX Va 



(15) 



where the right hand side is a constant that depends only 
on the grain size distribution. 

Thus, the quadron description makes it possible to col- 
lapse the statistics of all the systems, given the correct 
normalisation. As such, it gives better insight into the 
general characteristics of granular packs - characteristics 
that are independent of both the inter-granular friction 
coefficient and the initial state. In the next section, we 
explore the reasons for the apparent unified nature of the 
quadron volumes statistics. 



B. Cell order statistics 

To understand better the structural characteristics, we 
use a recently-proposed decomposition of the quadron 
volumes into conditional distributions [ll| 



P{v) = Y,Q{e)P{v 



(16) 



Here Q{e) is the occurrence probability of cells of order 
e and P{v\e) is the conditional PDF of the normalized 
quadron volume, given that it belongs to a cell of order e. 
Q{e) is essential to the understanding of random granular 
packing [l9| and it is this PDF that we wish to focus 
on next. As expected from Euler topological relation 
(below) and as can be observed from Figures |4]|6l a lower 
value of z must be accompanied with a higher mean value 
of e and corresponds to a lower packing fraction (j)' . 
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FIG. 12. The cell order probability Q(e) for for the 
10 systems generated from LIS and DIS. 

In Figure [12] we plot Q{e) for all the systems gen- 
erated from LIS and DIS. The plots make evident two 
points. One is that the higher the inter-granular friction 
the larger the fraction of high-order cells. The other is 
that (5(e) is hardly dependent on the initial state for any 
/i. 

There is a direct relation between the mean cell or- 
der, e, and the mean coordination number z and it can 
be derived from Euler's topological relation for a planar 
graph. 



Nv-Ne^Nc = 1 



(17) 



In this relation, Ny^ Ne and Nc are, respectively, the 
numbers of the graph's vertices, edges and cells. In this 
relation we disregard the one large cell making the out- 
side of the system, whose inclusion gives a topological 
characteristic 2, rather than 1, on the right hand side. 
This topological characteristic is negligible for sufficiently 
large systems. There are several ways to choose the graph 
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in granular assemblies. One is to take the grain centres as 
vertices and the lines connecting touching grains as edges. 
This construction is convenient for discs, but can be awk- 
ward for grains of arbitrary shapes. Another, more gener- 
ally applicable, construction is to take the contact points 
as vertices and the lines connecting contact points around 
grains as edges. This graph is consistent with the defi- 
nition of the quadrons (Figure [T]). Both graphs give the 
same result for packs of N{':^ 1) grains 
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where the rightmost term is a (negligible) boundary cor- 
rection. Figure [13] shows that the numerical systems sat- 
isfy this relation very well, showing that the sizes of the 
numerical pack are sufficiently large to ignore finite size 
effects. 
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FIG. 13. The mean cell order e vs. the mean coor- 
dination number z. 

Q{e) is sensitive to the value of ^, which is a function 
of /i (see Figure [12]), while P {v \ e) has been shown by 
Frenkel et al [11] (see also below) to be hardly dependent 
on /i. This suggests that there may be a parameterization 
of (5(e) that collapses all the curves corresponding to the 
collapse of P (u^) in Figure [TTl 

We find that all the curves collapse if we plot eeQ{e) as 
a function of e' = (e — e)/e^. We note that the integrals 
of the curves in Figure [M] are exactly 1. The collapsed 
curve appears to be fitted nicely by a truncated Gaussian 




eeQ{e) 



exp {-e'^ /2a^ll9) 



with erf{x) being the error function, e^^^ = 0.076 and 
(Jr = 0.082. However, we note small deviation from the 
Gaussian form at the large e-tail, for which we currently 
have no explanation. 



FIG. 14. The normalized cell probability eeQ{e). It 
is fitted well by a Gaussian form, except for small, 
but consistent, deviations in the large-e tail (high- 
lighted by a dash lined ellipse). 



C. Conditional quadron volume distributions 

To gain further insight into the universal properties 
of the structure, let us consider in more detail the con- 
ditional PDFs P{v\e). These PDFs were studied by 
Frenkel et al. [11], who argued, on the basis of geomet- 
rical considerations, that they should be independent of 
inter-granular friction. Their argument was based on the 
observation that, given a collection of N arbitrary grains, 
the number of ways to arrange e grains into a cell of or- 
der e depends only on the grains shapes and not on the 
inter-granular friction. While our results seem to provide 
a support to this argument, a closer look shows a system- 
atic /i-dependence of P {v \ e) for e > 6 (Figures [15] and 
[T6|) . We can also see this effect in the behaviour of the 
mean of the conditional quadron volume as a function of 
e, v{e) = ^q^^vP {v I e), shown in Figure [T71 We in- 
clude in the figure calculated values of quadron volumes 
of regular polygonal cells (RFC) of order e. 
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(20) 



which is clearly an upper bound for v{e). It is con- 
structive to consider the ratio of the observed mean 
quadron volume to the regular polygon value, 7(e) = 
Vq{e)/V^^^ (e), shown in Figure [T8l We observe that 
7(e), which is always below unity, has a minimum be- 
tween e = 5 and e = 6 and that it decreases with fi. 
These can be understood as follows. Carrying out a cell 
shape analysis (to be reported elsewhere), where cells 
were approximated to lowest order as ellipses, there is 
a strong correlation between 7(e) and the mean aspect 
ratio of the ellipses for any given e. The implication is 
that the increase of 7(e) with /i is due to the decrease 
of mechanical stability of elongated cells as /i decreases. 
The behaviour for any particular fi depends then on two 
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effects. As e increases, more cell configurations can be re- 
alised, allowing for more elongated cells and this accounts 
for the initial decrease of 7(e) with e. However, the more 
elongated the cell, the less stable it is, leading to fewer 
elongated cells as e increases. The competition between 
these two effects gives rise to the observed minimum. To 
test this understanding, we use next an analytical model. 




FIG. 15. The conditional PDF of the quadron vol- 
umes for e < 5 collapse nicely for all the systems. 




FIG. 16. The conditional PDF of the quadron vol- 
umes for e = 6, 8 collapse, but not as sharply as 
those for e < 5 (Figure [15)). 



V. A FREE CELL MODEL 

We believe that the geometry-based argument of 
Frenkel et al [11], mentioned above, should be augmented 
with stability considerations. Indeed, the number of ar- 
rangements of e grains, chosen from a given collection of 

grains, is independent of fi if we care only about the 
geometry. However, not all these shapes are stable un- 
der all possible combinations of compressive forces and it 
is this /i-dependent stability that constraints the possi- 
ble shapes. We demonstrate this effect, using a 'free cell 




4 6 
cell order e 

FIG. 17. The variation of the mean conditional 
quadron volume, with cell order e. The quadron 
volume of the regular polygon cell (dashed line) 
form an upper bound. 
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FIG. 18. The variation of the ratio of the mean 
quadron volume to regular polygon volume, 7 (< 
1), as a function of cell order e. 



model' (FCM). The model consists of a single cell, sur- 
rounded by e discs, chosen from the same distribution of 
sizes as in the numerical simulation. We construct a large 
number of cells for different values of \i and analyse their 
shape distributions. The cell construction is as follows. 
First, a random set of e discs is generated from the size 
distribution. We place the first disc at the origin and and 
another one touching it along the x-axis. We then gener- 
ate a random angle (^2 in Figure [19]) and place the third 
disc in contact with the second, making sure that it does 
not touch the first disc. This process is repeated, making 
sure that each added disc makes contact only with the 
previous one, until we reach the penultimate, (e — l)th, 
disc. We then determine the size and position of the last 
disc. If it fits in the remaining opening then the cell is 
kept, otherwise it is discarded and we attempt to gener- 
ate a new cell. Figure [20] shows examples of free cells, 
generated by this procedure. The FCM allows us to gen- 
erate quickly and efficiently a large number of cells and 
analyse them statistically. By comparing this geometry- 
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based model with the DEM results, we expect to identify 
the effects of mechanical stability for high order cells. 

A comparison between the resulting PDFs and those 
from DEM simulation for all e = 3 — 8 are shown in Fig- 
ure [2T1 The PDFs are similar, but with consistent differ- 
ences. The mean v for different e (Figure [T8|) are lower 
than of the numerical ones for e > 6. This systematic 
difference is the direct consequence of the cell stability 
constraint, which has not been considered by Frenkel et 
al. The DEM cells are closer to regular polygons because 
they must be stable under compressive forces. The more 
frictional the grains, the more elongated cells that they 
can support. This is the reason for the minimum in 7(e) 
decreasing and shifting to higher e with increasing fi. 




FIG. 19. The construction of a cell configuration 
studied by the free cell model. 




FIG. 20. Examples of free cell configurations for 
e = 6. 



VI. CONCLUSIONS AND DISCUSSIONS 

To conclude, we have studied in detail the effects of 
inter-granular friction and initial conditions on the struc- 
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FIG. 21. A comparison of the PDF of quadron 
volumes, generated both with the free cells model 
and with the DEM numerical experiments. 



tural characteristics of 2D granular assemblies. We have 
shown that structure can be analysed quantitatively, us- 
ing the quadron description [8, 9]. In particular, we have 
established that a number of distributions collapse onto 
single curves, pointing to a 'universal'-like structural 
characteristics that are independent not only of initial 
state, but also of the friction coefficient. These are sig- 
nificant results, whose implications are discussed at the 
end of this section. 

The raw quadron volume distribution, which is key for 
the statistical mechanics of granular matter [sl, [lol, HqI, 
changes systematically towards larger values with in- 
creasing friction coefficient fi and hence with decreasing 
mean coordination number z. However, when normalised 
by the mean quadron volume, it collapses nicely onto a 
single curve. Following the insight proposed by Frenkel 
et al.illl, we have traced this to a similar collapse of 
the conditional PDFs PiVq/Vg \ e) for every cell order 
e, except for small, but insightful, deviations for large e, 
which are discussed below. 

The common wisdom in the soil mechanics literature 
is that the relation between the the mean coordination 
number, and the packing fraction (j) depends on the 
initial state, as shown in Figure [3l However, prompted 
by the observation that the collapse appears once the 
rattlers are disregarded in the quadron partition of the 
granular space, we have examined this relation as a func- 
tion of rattlers- free packing fraction (j)' . We discover that, 
with this correction, the dependence on the initial states 
disappears and all the relations collapse onto a single 
curve, z{(j)'). This is consistent with our observation that 
the plot of Vq as a function of z is also independent of the 
initial state. Since rattlers do not support forces in our 
systems, these results carry a fundamental significance - 
they suggest that the packing fraction, the mean coor- 
dination number, the quadron volume distribution and 
generally the structure are linked directly to the manner 
in which stresses are transmitted. While this gives insight 
into the structure of the system, it should be remembered 
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that it is difficult to estimate the packing fraction of rat- 
tlers in real experiments and it is unclear at this stage 
how to use this insight for engineering applications. 

Considerations of mechanical stability also suggest 
that the argument proposed by Frenkel at el. [ll| , leading 
to independence of the conditional quadron volume dis- 
tributions of inter-granular friction, needs to be modified. 
Their argument assumes that cell shapes are independent 
of mechanics and therefore that arranging an e-sided cell 
depends only on the grain shape distribution. This, in 
turn, eliminates friction as a factor in the occurence prob- 
ability of cell structures. However, as friction decreases, 
elongated cells are less stable mechanically, reducing the 
number of available cell shapes. Indeed, we can observe 
a small increase in the occurrence probability of large 
quadron volumes in P{Vq/Vg \ e) for e > 6 as /i increases 
(Figure [T6|). 

This effect is seen most clearly when plotting the ratio 
of the mean quadron volume to that of a regular polygon 
as a function of e (Figure [T8|). The initial drop reflects 
the departure of the cell from a regular polygon due to 
the increased number of possible geometrical configura- 
tions, while the subsequent increase shows the effect of 
mechanical stability in limiting such a departure. The 
more pronounced increase for lower /i is evidence for our 
argument above. 

To illustrate the effect of mechanical stability, we have 
constructed a free cell model (FCM) that takes into ac- 
count only geometric considerations and disregards me- 
chanical stability. Indeed, this model shows the same 
initial decrease of the above ratio with e and none of 
the subsequent increase (Figure [18]). Thus, neglecting 
the mechanical stability effect, the FCM provides a low 
bound for the mean quadron volume, as illustrated in 
Figure [221 
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FIG. 22. Simple prediction of the mean quadron 
volume in terms of the mean coordination number. 

We have further found that the raw cell order PDF, 
(5(e), changes systematically with /i, but hardly at all 
with the initial state. However, all the PDFs also col- 
lapse to a single curve under the transformations Q{e) 
eeQ{e) and e (e — e) /e^. We note that the collapse 



cannot be perfect at the large e tail due to the aforemen- 
tioned constraint of mechanical stability. This friction- 
specific behaviour, which is evidenced by Figures [16] and 
[T8l is hardly noticed at all in Figure [HI demonstrating 
that the effect of mechanical stability is small compared 
to the geometrical effect discussed by Frenkel et al. |Tl|. 

Our results are significant for several reasons. First, 
they are a step towards a systematic understanding of 
structural organisation of granular matter in response to 
a specific packing procedure. In particular, they suggest 
that the structural characteristics of granular matter are 
best understood via the statistics of cell configurations. 
Second, they make it possible to progress on the statisti- 
cal mechanics of granular matter, where knowledge of the 
quadron volumes distribution is key. Third, the results 
shed light on the common wisdom in the soil mechanics 
and civil engineering communities that the initial state 
affects the final relation between z and ^. 

Our findings show that this relation is directly linked 
to mechanical stability of the packing backbone - for 
any specific pack generation process, the initial state af- 
fects the fraction of rattlers, which do not participate 
in the stress transmission. Once the rattlers are disre- 
garded, there emerges a unique relation between z and 
the corrected rattlers- free packing fraction 0^ Fourth, 
our analysis has a significant ramification for the famous 
packing problem. Different packing procedures lead to 
different structures. For example, procedures not con- 
strained by mechanical stability would give rise to dif- 
ferent (5(e) 's than those that are. Specifically, the for- 
mer would tend to have lower packing fractions due to 
a higher occurrence probability of large cells, which are 
mechanically unstable. This reinforces the need for a 
more accurate definition of the packing problem, which 
we discuss elsewhere [19]. 

Taken together, these results show that the pack- 
ing process has a fundamental effect on the structural 
characteristics while the initial conditions and the inter- 
granular friction are details that only modify this effect 
and can be scaled away. The collapsed distributions show 
how to predict such modifications. It should be empha- 
sised that the collapses, which we have found, do not im- 
ply universality in the traditional sense because the spe- 
cific forms of the collapsed curves depend on the packing 
procedure. For example, it has been demonstrated that, 
using a different packing protocol, it is possible to gener- 
ate packs dominated by cells of order 3 [20], which can- 
not be scaled to match the ones we obtained here. There 
is in principle an infinite number of possible procedures 
to generate packs, depending on an astronomically large 
number of parameters. It would be impossible to map all 
the possible processes to a manageable parameter space. 
Each such procedure gives rise to its own characteristic 
(5(e) and it would be impossible to collapse all of these 
(5(e) onto one curve. 

To fully understand how structural characteristics de- 
pend on the packing procedure, one needs a good model 
of both the geometric effects and the limitations that 
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mechanical stability poses on the cell shape distribution. 
We are currently developing such a model and will report 
results in a later publication. Intriguingly, the relevance 
of both volumetric and stress effects is reminiscent of the 



inter-dependence demonstrated recently for the statisti- 
cal mechanical understanding of granular matter ^| and 
it would be interesting to find out whether the two issues 
are related or not. 
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